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Abstract 

Several hadronic observables have been studied in central 158^ GeV Pb+Pb collisions 
using data measured by the WA98 experiment at CERN: single vr" and K~ production, as 
well as two- and three-pion interferometry. The Wiedemann-Heinz hydrodynamical model 
has been fitted to the pion spectrum, giving an estimate of the temperature and transverse 
flow velocity. Bose-Einstein correlations between two identified vr" have been analysed as a 
function of kx, using two different parameterizations. The results indicate that the source 
does not have a strictly boost invariant expansion or spend time in a long-lived intermediate 
phase. A comparison between data and a hydrodynamical based simulation shows very good 
agreement for the radii parameters as a function of k^- The pion phase-space density at 
freeze-out has been measured and agrees well with the Tomasik-Heinz model. A large pion 
chemical potential close to the condensation limit of 771,,^ seems to be excluded. The three- 
pion Bose-Einstein interferometry shows a substantial contribution of the genuine three-pion 
correlation, but not quite as large as expected for a fully chaotic and symmetric source. 

1 Introduction 

The study of single particle distributions of particles produced in heavy-ion collisions gives 
access to the degree of thermal and chemical equilibrium at freeze-out and allows the deter- 
mination of the parameters of hydrodynamical expansion models of the source. 

The spatio-temporal extension of the interaction region created in such collisions is not 
directly observable, but the study of Bose-Einstein interferometry between identical particles 
provides information on the geometry and on the dynamical evolution of the particle emission 
sources. In particular, the correlations between produced pions give access to the size of the 
homogeneity region, to the duration of emission, and to various parameters characterizing 
the spatial extension of the fireball Q. In addition, by combining information from the source 
size in momentum-space obtained by interferometry, and from the momentum-space density 
provided by the single particle distributions, an average phase-space density at freeze-out 
can be calculated. 

Compared to the two-particle correlation, the three-particle correlation can provide ad- 
ditional information on the chaoticity and asymmetry of the source emission ||2|, |^, In 
particular, the three-pion interference produced by a fully chaotic source is sensitive to the 
phase of the Fourier transform of the source emission function and, hence, to the asymmetry 
of the source. 

In this paper, we present the analysis of single particle production, and of two- and three- 
pion interferometry measured in the WA98 experiment for central 158A GeV 208p|j_j_208p)2 
collisions at the CERN SPS. In addition, estimates of hydrodynamical expansion model 
parameters, temperature and transverse flow velocity are extracted from these results. 

2 Experimental setup and data processing 

The CERN SPS fixed target experiment WA98 p combined large acceptance photon de- 
tectors with a two arm charged particle tracking spectrometer. The experimental layout is 
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shown in Fig. |^. The 158A GeV Pb beam interacted with a Pb target near the entrance 
of a large dipole magnet. The onhne trigger centrahty selection used a forward calorimeter 
located at zero degrees and a mid-rapidity calorimeter measuring the total transverse energy 
in the pseudorapidity interval 3.2 < rj < 5.4. The results presented here have been obtained 
from an analysis of the complete data set. These data were taken with the most central trig- 
gers corresponding to about 10% of the minimum bias cross section of 6300 mb |^], with 
an average of 330 participating nucleons per collision. These quantities are estimated to have 
an overall systematic error of less than 10%. The charged particle spectrometer made use 




Figure 1: The WA98 experimental setup. 



of a 1.6 Tm dipole magnet with a 2.4x1.6 m^ air gap which deflected the charged particles 
in the horizontal plane into two tracking arms located downstream, one on each side of the 
beam axis. The results shown here are measurements of vr" and K~ observed in the nega- 
tive particle tracking arm of the spectrometer. This tracking arm consisted of six multistep 
avalanche chambers with optical readout Inside the chambers, triethylamine (TEA) pho- 
toemissive vapour produced UV photons along the path of charged particles, these photons 
being subsequently converted to visible light via wavelength shifter plates. On exit, the light 
was reflected by mirrors at 45° to CCD cameras equipped with two image intensifiers. The 
active area of the first chamber was 1.2x0.8 m^ and that of the other five 1.6x1.2 m^. Each 
CCD camera pixel viewed a region of about 3.1x3.1 mm^ on the chambers. Downstream 
of the chambers, at a distance of 16.5 m from the target, a 4x1.9 m^ time of flight wall 
allowed for particle identification with a time resolution of better than 120 ps. The resulting 
particle separation is shown in Fig.^. The tt~ rapidity acceptance ranged from y = 2.1 to 
3.1 with a rapidity average at 2.7, close to the mid-rapidity value of 2.9. The momentum 
resolution of the spectrometer was Ap/p = 0.005 at p = 1.5 GeV/c, resulting in an average 
precision of better than or equal to 10 MeV/c at vanishing px for all the Q variables used in 
the correlation analysis and defined in sections ^ and ^: (7{Qinv) = 7 MeV/c, (t{Qto) = 10 
MeV/c, a{QTs) = 5 MeV/c, a{QL) = 3 MeV/c, c7(Qt) = 8 MeV/c, a{Qo) = 5 MeV/c, 
a{Q^) = 7 MeV/c. 

Severe track quality cuts were applied, resulting in a final sample of 7.9 x 10^ vr^, providing 
13.7x10^ pairs and 13.1x10^ triplets for the Bose-Einstein correlation analysis. 
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Figure 2: Particle identification in the negative arm of tlie spectrometer. 



3 Single particle spectra and hydrodynamical model 

The study of inclusive distributions of single particles produced in heavy-ion collisions can be 
interpreted in the context of models of the source using hydrodynamical expansion. Within 
the context of such models, which assume local thermal equilibration, parameters like the 
temperature and collective velocity at freeze-out can be determined. In the limit of a sta- 
tionary fireball, the distribution takes the simple form Q 

E = oc Ee~^ (1) 

dp^ mTduiTdydcj) 



where p is the Cartesian particle momentum, mx = yPj- + "^o transverse mass, pT is 

the transverse momentum, mo is the rest mass, y is the rapidity, E = mTCOsh.{y — y fireball) , 
jjL is the chemical potential and T is the temperature. In the limit where only a narrow 
rapidity interval close to y fireball is measured, the spectrum becomes 

dp'' 

and in the case of a rapidity-integrated spectrum 

mTdrriT \ T J 

where Ki is a modified Bessel function. The last approximation holds in the limit mj- ^ 
T. Plotted against mx — mo, all particles from a thermalized emitter should show the 
same universal exponential behaviour. However, different additional features like transverse 
hydrodynamic expansion or particles originating from the decay of resonances will distort the 
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shape of the spectrum. It is noted in |1C] that for the popular fit with a simple exponential 
in rriT (without the mT-prefactor) 

= Ce-^'T/T (2) 

an interpretation of the resulting slope parameter in terms of a temperature is not possible. 
However, since it is found to fit the measured spectrum better than the previous expressions, 
it is useful for obtaining an estimate of the inverse slope parameter T. 

3.1 Data analysis 

A detailed description of the analysis of the single particle spectra presented here can be 



found in |11]. The correction for detector acceptance and efficiency applied to the measured 
spectrum is based on a precise simulation of the detector, tuned in order to reproduce the 
measured detector response as accurately as possible using the VENUS [^] event generator 
as input. Multiple scattering, decays and all other reactions within the detector material 
are taken into account. Efficiency maps depending on the hit position on each chamber, the 
particle momentum and its identity are applied, together with position resolution and noise 
simulation. Simulated events are then reconstructed using the same code as for real data, 
ensuring that any software- induced systematics are also taken into account. The output is 
then matched to the VENUS input. This correction procedure ensures that the final result is 
little sensitive to remaining contamination, such as pions in a kaon sample. A detailed study 
of the systematic uncertainty has been performed and it is found to be small, especially as 
regards the slope of the spectrum (~3%). The absolute normalisation on the other hand is 
more sensitive to detector instabilities which could not be perfectly simulated and is found 
to have an estimated uncertainty of at most 20%. 

Since the statistical uncertainty on the final mj- spectrum is negligible compared to the 
systematics, it is possible to apply very severe quality cuts on the reconstructed tracks and 
on the events used. In contrast to the Bose-Einstein analysis, only a subset of all data is 
kept using only run periods where the detector operation remained particularly stable. The 
analysis is performed separately for identified 7r~ and K~. The final data sample consists 
of 4.7x10^ 7r~ tracks and 1.8x10^ K~ drawn from 3.8x10^ central events. Figures ^ and ^ 
show the detector acceptance for negative pions and kaons respectively, obtained from the 
simulation. The Tr^/vr^ ratio and the K~ /K'^ ratio of the number of detected pions and 
kaons have also been measured using two opposite magnetic field polarities. The results will 
be presented in a separate publication. 

3.2 Results 

Figure ^ shows the final (fully corrected) single particle spectrum for 7r~ plotted as a function 
of the two variables mx and y. The form of the detector acceptance shown in Fig. ^ is 
clearly noticeable, except for the low py-low y lobe, which has been omitted here. Figure 
^ shows the projection on the mT axis (normalized to a unit rapidity interval using the 
VENUS profile) of the two-dimensional spectrum of Fig. A fit to the exponential form 
of eq. |2| over the interval rriT — ttiq = [0.1,1.2] GeV/c^ (not shown in Fig. |) yields C = 
5220 ± 80 {stat.) t^g^" {syst.), and T = 0.168 ± 0.001 [stat.) tg;[Jg| [syst.) GeV, with a 
per degree of freedom of 1.6. This result is in agreement with the NA49 result [13|. 
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Figure 3: Acceptance of tracking arm I in the {pT,y) plane for pions. 
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Figure 4: Acceptance of tracking arm I in the (pr, y) plane for kaons. 
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Figure 5: Two-dimensional (fully corrected) single pion spectrum as a function of the two variables 
rriT and y. 



At high pT, where perturbative QCD becomes applicable, the spectra are expected to 
attain a power-law behaviour as observed in numerous high energy pp measurements (see 
for example The heavy- ion data of this experiment seem to follow that trend even 

into the lower px range. Therefore, a parametrization originally inspired by QCD ||l5| and 



successfully applied already to pp data [14| and heavy- ion data |16] can be used to fit the 
spectrum: 



dp-^ \PT+PoJ 

with C, po) and n taken as free parameters. A link to the more familiar exponential slope 
parameter T is obtained from the derivative of this expression according to 



f{PT) _Po , PT 

dpT 



J power-law - " " " + ^ 



Thus, p^jn characterizes the slope of the transverse momentum spectrum in the limit pT 0, 
while 1/n characterizes its gradient along p^, i.e. the strength of the concave curvature. The 
extracted parameters are C = 4150 it 70, po = 4.80 it 1.04 GeV/c and n = 29.0 it 5.9, which 
gives a slope parameter T = po/n = 0.166 ± 0.005 GeV. The P^r degree of freedom is 1.0. 
The same fit performed on the vr^ spectrum measured by the WA98 experiment ||^, ^, 17, 18 1 



yields C = 5120±140, po = 5.08±0.18 GeV/c, n = 29.3±0.8, and so T = 0.173±0.002 GeV. 
The acceptance for the measurement being different, it is interesting to note that fitting 
only the pT = [0.3, 1.4] GeV /c common interval of both spectra yields T = 0.169 GeV for 
the 7r~ and T = 0.166 GeV for the vr*^, showing that the result is stable with respect to the 
fit interval and that the spectra agree quite well. 
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It has also been shown that fluctuation of the parameter in an exponential distribution 
leads to a final distribution of the power-like form [^]. This behaviour can be interpreted 



in terms of a suitable application of the nonextensive statistics of Tsallis [20]. This interpre- 
tation is convenient to describe particle production at fluctuating T, as may occur near the 
phase transition. These fluctuations would exist in small parts of the hadronic system with 
respect to the whole system rather than between events. The average (T) around which the 
temperature fluctuates is given hy po/n = 0.166 GeV and the relative variance of 1/T is 



- = 0.034 , 

n 



both for TT~ and vr''. This corresponds to a nonextensivity parameter q = l-\-uj of 1.034. This 
result, interpreted in the spirit of |]l9|, indicates a relative fluctuation AT/T of 18.4 it 1.9% 
(18.4 lb 0.3%) for the vr" (vr'^) measurement. However, this analysis neglects the contributions 
of resonance decays and of flow velocity distributions. Both of these effects increase the 
curvature of the pion spectrum. Thus the value of AT/T = 18.4% should be considered an 
upper limit on the temperature fluctuations. 

The averaged negative pion yield per unit rapidity in the acceptance window is 1/NcvtdN/dy 
129 lb 1 (stat.) tl^ (syst.). Figure shows the same plot as above for kaons. A fit 
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Figure 6: One-dimensional (fully corrected) mj-spectrum for tt . The errors are statistical. The 
power-law fit of eq. ^ is superimposed on the data points. The fit interval is rriT — tuq = 
[0.05, 1.2] GeV/c^. 



to the exponential form | yields C = 1000 ib 120 (stat.) If (syst.), and T = 0.181 ib 
0.005 (stat.) iE):oo9 (syst.) GeV. So the inverse slope T for kaons and pious are comparable. 
It should be noted that the rapidity acceptance is quite different for the two particle species, 
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Figure 7: tt^ (black circles) and K~ (empty circles) single particle spectra. The error bars are 
statistical only. Note that for the same transverse mass, the pt is different for pions and kaons. 
The inset shows the (bin by bin) ratio of the two px spectra. 



the TT~ acceptance being much closer to mid-rapidity. An estimate of the effect induced 
by this difference in acceptance can be made using the term cosh(Ay) given below eq. |l[ 
The inverse slope for K~ would become of the order of 0.230 GeV. The /tt^ ratio at a 
common rapidity of 2.3 is 11.4 it 0.4 (stat.) t":? {syst.)%, which is in agreement with the 
NA49 result Q. 



3.3 Hydrodynamical Source Expansion Model 

The single pion spectrum has been fitted in the interval 0.05 < niT — mo < 1.2 GeV/c^ with 
the Wiedemann-Heinz (W.-H.) model ||2l|, 23, 24, 25, 26|. It relies on the following idea: 
the main characteristics of the particle phase-space distribution at freeze-out can be quanti- 
fied by its widths in the spatial and temporal directions, a collective dynamical component 
(parameterized by a flow field) which determines the strength of the position-momentum 
correlations in the source, and a second, random dynamical component in momentum space 
(parameterized by a temperature). The model emission function contains seven parameters, 
but the shape of the single particle transverse mass spectrum is fully determined by the tem- 
perature T and the transverse flow rapidity profile ryr(r) = Vf^ which is assumed to depend 
linearly on the transverse coordinate r, where f// is the transverse flow rapidity strength, 
and R the Gaussian transverse spatial width. The mean transverse flow velocity (/3t) can 
be easily calculated as the mean value of tanh(r7x'(r)) over the transverse source profile. The 
results will be given as a function of rather than r]f, since its physical interpretation is 
more straightforward. 

The three parameters that are allowed to vary freely during the fitting procedure are T, 
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r]f, and an overall normalization factor. The result of the fit is T = 0.084 it 0.003 GeV and 
(Pt) = 0.50 ± 0.02, the x^/d-o.f. being 1.1. The resulting curve is not shown in Fig. |6| as 
it is hardly distinguishable from the power-law fit. The temperature and flow parameters 
are strongly correlated, so it is more interesting to consider a contour plot of the fit 
to the measured single particle spectrum as a function of those two parameters. Figure ^ 
shows the result assuming a Gaussian transverse spatial profile, and Fig. ^ assuming a box 
profile of the same rms width as the Gaussian one. Only direct pions are considered in the 
model. Furthermore, only statistical errors are taken into account during the fit. The curves 
displayed represent different confidence levels for the (T, (Pt)) values. The dependence of 




0.065 U ... I ... I ... I ... I ... I .. I 
0.44 0.46 0.48 0.5 0.52 0.54 

(Pt> 

Figure 8: contour plot of the hydrodynamical model fit to the measured vr" single particle 
spectrum of Fig. |^. A Gaussian profile is assumed for the transverse spatial profile of the source. 
The curves displayed represent successively (starting with the innermost one) contours at 50% 
CL, 90% CL, 99% CL (dashed line), 99.9% CL, and the last curve indicates a highly excluded 
region. The black square represents the best parameter values. 

the x^/d-o.f. on the shape of the source is very small, the best x^/d-o.f. being slightly larger 
for the box profile (1.5). While the best fit temperature appears to be independent of the 
shape of the source, the best fit fiow velocity is considerably smaller ((/3t) = 0.29) for the 
box profile. Inclusion of the resonance decay contribution in the model calculation would 
be expected to reduce the extracted flow velocity parameters and increase the temperature 
parameters. The precise effects do however depend on details of the models used. 

4 Two-pion correlations 

Bose-Einstein interferometry is most commonly used to study pairs of identical particles. 
The correlation function C2 is deflned, up to a proportionality factor M, as the ratio of a 
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Figure 9: Same plot as in Fig. ^ but assuming a box profile for the transverse spatial profile 
the source. 



two-particle spectrum V2 over the product of two single particle spectra Vi: 



^2(Pl,P2) 

Vi{pi) ■ Vl{p2 



with 



and 



Vi{p)=E 



dN 
dp^ 



'P2{pi,P2) = El E2 



dN 



dPi dp2 

where Ei and pi are the energy and momentum of particle i, respectively. 

Experimentally, the product of one-particle distributions in the denominator is commonly 
obtained by a mixed event technique whereas the two-particle distribution in the numerator 
is constructed from all pair combinations of identical particles found in each event. C2 is 
then normalized to unity far away from the interference region. 

A fully chaotic source can be seen as a superposition of uncorrelated elementary sources, 
and one- and two-particle distributions may be expressed through the Wigner function of 
the source S{x,ki2) Q. The correlation function is then written 



C2{Pl,P2] 



1 + A 



/ d'^x S{x,ki2) exp[igi2x] 
\!d^x 5(x,fci2)|2 



with qi2 = pi — P2, the 4-momentum difference of the two particles, ku = (pi + p2)/2 and 
X = (xi + X2)/2, the mean space-time coordinate of the pair emission point. The chaoticity 
parameter A is inserted to take into account the possibility that the source may not be fully 
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chaotic and also that any wrongly reconstructed tracks, or tracks originating from decays of 
long-lived resonances, will dilute the Bose-Einstein correlations in the data. 

A one-dimensional interferometry analysis is commonly made as a function of Qmv = 
\J—qi2, whereas a multidimensional analysis uses a set of Q variables which are defined as 
various projections of qi2- 



4.1 Data analysis and results 

Two independent analyses were performed with the complete data set recorded in the nega- 
tive tracking arm [0, totalling 13. 7x 10^ pairs of identified tt~ . These data were corrected 
for resolution and Coulomb effects in an iterative way [p8| . The Gamow correction was not 
used as it overcorrects the data in the Qinv range of 0.1 to 0.3 GeV/c. Because the finite 
resolution in the measurement of the Q variables leads to an underestimate of the radii and 
A parameters, a correction has to be implemented in the fitting procedure. This is done by 
a convolution method, replacing the C2{Q) formula expressing the two-particle correlation 
function used to fit the data by 

C^2'(Q) = /// riQ,Q') dQ' 

where r{Q,Q') is the resolution function which is chosen to be Gaussian: 
riQ,Q') = l/i27rf/^ l/|y|i/2 

X exp[-l/2 (Q - Q'f (Q - Q')] 

The diagonal elements of the covariance matrix V are set equal to the square of the various 
Q-resolutions, which are estimated by a full simulation of the experimental setup as a function 
of the kr = \pti +PT2I/2 of the pairs. The non-diagonal elements of V are neglected, as the 
resolution correction has a very small effect compared to the Coulomb correction. Figure 



10 shows C2, the measured 7r~7T~ correlation function, plotted as a function of Qinv before 

and after the resolution and Coulomb corrections. In the plot, the correction for resolution 

is obtained by multiplying each data point by C2^{Qinv) / C2{Qinv) ■ C2 is clearly exponential 

p9| . The solid curve is a fit of the form l+Xeexp[-2QinvRinv] which gives Rinv = 7.33 ib 0.08 

fm and Ag = 0.788ib0.009 for (A;t) = 0.116 GeV/c. This exponential behaviour appears not to 

hold in the 3-d analysis, where the projected slices are better represented by Gaussians. The 

3-d analysis of Bose-Einstein correlations has been done using two different parameterizations 

in the Longitudinally CoMoving System (LCMS): the standard Pratt-Bertsch fit (PB) in the 

3-dimensional space of momentum differences Qts (perpendicular to the beam axis and to 

the transverse momentum of the pair), Qto (perpendicular to the beam axis and parallel to 

the transverse momentum of the pair), and Ql (parallel to the beam axis) pO|, including a 
2 

out— long 



cross term /„„„ as predicted 



C2 — 1 + Xexp[—Q'^gR^g — Q^qR^q — Q\R\ 

-2QtoQ iRlut-long] 

and the Yano-Koonin-Podgoretskii fit (YKP) |^l| in the Qo (energy difference of the pair), 
QtiQl space 

C2 = 1 + Aexp[-Q^i?f. + {Ql - Ql)Rl - {Q-Uf{Rl + i?|)] 
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Figure 10: The measured two-pion correlation function C2 (full symbols), corrected for resolution 
and Coulomb effects, as a function of Qmv The full line is a fit to an exponential form whereas 
the dashed curve is a fit to a Gaussian form. The empty symbols show the data before corrections. 
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Figure 11: The kx dependence of the Pratt-Bertsch parameters. The points are plotted at the 
average kx of the bin and the horizontal bars indicate the bin width. The open symbols in the 
figure at bottom right show the A parameter after correction for background from misidentified 
pions, as explained in section ^. 
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where U = 7(1,0, 0,Ui), 7 = — vf^ with vl in units of c = 1. In the YKP parame- 

terization the different radii are invariant under a longitudinal Lorentz boost, and the speed 
parameter vl connects the local rest frame to the measurement frame (the LCMS in our 
The results as a function of the of the pairs are shown in Figs. O and 12 and 



case 



summarized in Table The A parameters of the YKP fit (not shown in Table |^) are found 
compatible with the A parameters of the PB fit. 
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Figure 12: Comparison between the PB and YKP fits in the LCMS. The cross term i?; 
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out- 
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from 



the PB fit and vl from the YKP fit both show a deviation from zero at large kx when estimated in 
the LCMS. The Rts and Rl parameters from the PB fit are in good agreement with respectively 
Rt and R^ from the YKP fit. 



The systematic errors, not included in Figs. 11 and 12, and Table |l], are estimated by 
varying the different analysis cuts, including the cuts used to identify the pion with the time 
of flight system. The systematic error on the Coulomb correction due to the error on the 
determination of the radius parameters is also taken into account. All these variations are 
added in quadrature. The total relative systematic errors on the radii AR/R amount to 
0.8% for Rinv, 1.4% for Rts, 3.5% for i?TO, 9.1% for Rl, 0.8% for Rt, and 9.7% for R^. 
Systematic errors on R'^uf—iong '^nd V£, are asymmetric and reach respectively i?;? fm^ and 



+0.13 
-0.08 



The Rts and Rl parameters from the PB fit are in good agreement respectively with 
Rt and i?4 from the YKP fit. The cross term R^ut-iong froni the PB fit and vl from the 
YKP fit deviate from zero. In a source undergoing a boost invariant expansion the local rest 
frame coincides with the LCMS. Both the cross term and vl expressed in the LCMS are then 



expected to vanish |21|. As this is not quite the case, it suggests that the source seen within 
the acceptance does not have a strictly boost invariant expansion. The strong decrease of the 
longitudinal radius Rl or R4 with kT compared to the transverse radii Rt, Rts, Rto shows 
a longitudinal expansion which is larger than the transverse one. Finally, the Rq parameter 
(not shown in the figures), which corresponds to the duration of emission of particles from 
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the source, is compatible with zero for ah bins, excluding a long-lived intermediate phase. 
These results agree with the previously published WA98 results obtained using roughly half 
of the present data sample p9|. The WA98 analysis is compatible with the NA49 results 
obtained in a slightly different rapidity range ((y) =3.2) with unidentified negative particles 



4.2 Comparison with a hydrodynamical model based simula- 
tion 

The W.-H. model can be used to generate correlation functions which can be compared to 



the data. For that purpose, all necessary integrations are performed numerically [25| to get 
the value of C2 in the PB parameterization for given values of Qts^ QtO: Ql, ^Tj and the 
rapidity Y = {yi + 1/2)/'^ of a pair. To simulate properly the acceptance of the negative 
tracking arm, the mean values of kx and Y for each {QtSiQto,Ql) bin are calculated 
using real data, and then for each bin the mean values {{Qts) j {Qto) AQl), {kr) , (Y)) are 
used as input for the hydrodynamical calculation. This procedure is repeated for all five 
kx intervals used in the data analysis. The correlation functions are generated neglecting 
contribution from resonances, using T = 85 MeV, r]f = 0.5 (values extracted from the single 
particle spectra analysis), R = 8 fm, tq = 11 fm/c, Ar = 2 fm/c and Arj = 1.3 (see [ ^ ] 
for definitions of the parameters). Figure |l^ shows the Bose-Einstein radii extracted from 
fitting the simulated correlation functions with the PB formula. The error bars used to 
perform these fits are taken to be the same as the ones calculated for the real data in each 
(Qts, Qto, Ql) bin. The measured results are shown on the same plot, and the agreement 
is found to be good. The shift in the the cross term Rgut-iong between data and simulation 
is compatible with the systematic uncertainty. 



5 Average pion phase-space density at freeze-out 

As the mT spectrum gives the momentum-space density at freeze-out and as the Bose- 
Einstein correlation radii provide information on the covariant volume for particles of mo- 
mentum p, it is possible, by combining these results, to extract the average phase-space 
density (/) {pT,y) at freeze-out [^, |3^ : 



"v/A dn 



with = \J + = my cosh y. The radii and A are functions of pT and y. The factor 
A, which comes from the two-pion correlation analysis, corrects for contributions of pions 
originating from long-lived resonances decaying close to the primary vertex. A difficulty 
of this method is to include only the contribution of the real pions in the determination 
of A, excluding backgrounds from misidentified particles. This is achieved by applying to 
A, separately for each kx bin, a correction factor obtained from a full simulation of the 
experimental setup, taking into account geometrical acceptance, backgrounds and efficiency 
of the chamber-camera-time of flight system. The effects can be seen in Fig. O, bottom 



right, where A with and without correction is displayed. Figure |lj and Table y show the 
results on the average phase-space density for vr" as a function of pT- The error bars reflect 
the statistical errors only. The systematic uncertainties are dominated by the uncertainty 
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[GeV/c] 



Figure 13: Comparison between simulation and data as a function of kx- The open symbols are 
the result of the fit of the PB formula to the correlation functions produced by the hydrodynamical 
model, whereas the full symbols are the result of the fit to the data. 
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Figure 14: Average pion phase-space density as a function of pr compared to different models (see 
text). 
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of the correction on the A parameter, which is estimated to be 20%, and by the systematic 
error on i?/,, giving a total of 13.7% systematic error on (/). Within errors, all nuclear 
collision measurements at the SPS are found to be indistinguishable |35|, and our result at 
mid-rapidity agrees well with these previous measurements. The dashed lines indicate Bose- 
Einstein density distributions (/) = [eyiY>{E-n /T) — 1]^-*^ of static sources of pions {E.,^ w my) 
for three choices of the freeze-out temperature T: 80 MeV, 100 MeV, and 120 MeV, from the 
lower to the higher curve. The results are in rough agreement with the 100 MeV distribution 
at low pt, but show a clear deviation from a Bose-Einstein distribution at high pT- As pointed 
out in [^6|| , this deviation is mainly due to the strong longitudinal expansion of the fireball 
which reduces the spatially averaged phase-space density and, to a lesser extent, due to the 
radial collective flow, which adds extra transverse momentum to the particles compared to 
particles emitted by a static source. Consequently, even in the absence of transverse flow, (/) 
will be reduced compared to a Bose-Einstein density distribution. This effect may necessitate 
a positive pion chemical potential in order to match the experimental observation. Such a 
positive potential can be related to the presence of pions from short-lived resonance decays. 
The Tomasik-Heinz (T.-H.) model |3^] used to fit the data assumes a thermalized fireball 
with a longitudinally boost-invariant expansion and a transverse fiow rapidity profile which 
depends linearly on the transverse coordinate. This model includes a pion chemical potential 
and has three free parameters, the freeze-out temperature T, the strength of the transverse 
flow rapidity profile r/j, and /io, the chemical potential value in the center of the fireball. In 
Fig. 14, the full line is a fit of the T.-H. model with a box transverse density profile, whereas 
the point-dashed line is a fit of the same model with a Gaussian transverse profile. Both fits 
agree well with the data with a /d.o.f. of 0.78 and 0.93 for respectively the box and the 
Gaussian proflles, giving T = 8711? MeV, 77* = 0.49 tg:^^ and ^uq = 571?^ MeV for the 
box profile. This r]t result corresponds to a mean transverse fiow velocity {(St) = 0.42±o.il • 
Due to the lack of experimental points at large pTt the T.-H. model with the Gaussian profile 
provides very loose estimates of these parameters, which are nevertheless compatible with 
those from the box profile. Figure 15 shows the contour plot for the T.-H. model with 
the box profile for T and r/j. The curves represent (starting from the centre) contours at 
39%, 70% and 99% confidence level. The best fit is obtained for a pion chemical potential 
/xo = 57 MeV (black square) but a solution with /xq = MeV is also possible (full circle). On 
the other hand, a large pion chemical potential at freeze-out, close to the Bose condensation 
limit of iiQ = m^, as could be speculated given the rather low mass of the pion, seems to 
be excluded in view of the error bar on /io- The results of the fit of the W.-H. model on 
the single vr" spectrum can be compared to the fit of the T.-H. model on the phase-space 
distribution. The agreement is good for T, and satisfactory for {(3t) when taking into account 
the systematic error. It should be noted that the constraint provided by the fit of the T.-H. 
model is weak compared to the one given by the W.-H. model because of the relatively small 
amount of experimental points in the phase-space distribution. Moreover the T.-H. model 
uses a pion chemical potential whereas the W.-H. model doesn't. Finally, as (/) is the pion 
occupation per 6-d position(8)momentum cell, the obtained values do not provide striking 
evidence for the presence of an excess of pions or for the presence of large disoriented chiral 
condensates. This chiral condensate phenomena has been investigated by other means within 
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Figure 15: contour plot for the T.-H. model with a box profile. Starting from the centre, the 
curves represent contours at 39%, 70% and 99% confidence level for T and r\t. The black square 
is the best fit (//q = 57 MeV). The full circle is the result with /jq excluded from the fit and set to 
zero. 
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6 Three-pion correlations 



In the hypothesis of a fully chaotic source of identical particles, the two-pion correlation 
function can be written C2 = 1 + |-^i2p where I-F12P is the Fourier transform squared of the 
space-time source function. The three-pion correlation function is then C3 = 1 + |i^i2p + 
I i^23 P + I i^3i P + 2 • Re{Fi2 • F23 • F31 } . The terms | Fij \ ^ , which express the contribution of the 
three combinations of the two-pion correlations contained in the triplet (123), provide the 
largest contribution to C3. The last term only represents the genuine three-body correlation. 
It can be written 2 • |Fi2| • IF23I • IF31I • W, where Fij is defined as \Fij\ ex.p[icj}ij] and W = 
cos((/)i2 + </'23 + 031 )• The simultaneous measurement of C2 and C3 provides information on 
W, the cosine of the sum of the three phases of the Fourier transforms. In contrast, the 
measurement of C2 alone gives access only to the radii of the source and not to the phases. If 
the emission source is fully chaotic, W measures the asymmetry of the source. In the presence 
of not fully chaotic sources, which is likely to be the case, W , the strength of the true three- 
body correlation, is basically sensitive to the coherence. So, in the framework of the partially 
coherent model [^] , W gives information on the degree of chaoticity of the emission source 
in a manner which is insensitive to backgrounds such as the contributions from resonances. 
The complete data set recorded with the negative tracking arm yielded a total of 13.1x10^ 
triplets of vr". After correction for resolution and Coulomb effects Q a strong C3 signal is 
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Q3 [GeV/c] 

Figure 16: The three-pion correlation function C3 as a function of Q3. The full line is a fit to 
a double exponential form (see text). The dotted and the dashed lines are fits to respectively a 
single exponential and a Gaussian form. 



observed (Fig. ^) as a function of Q3 = yQi2 + Q23 + Q31 with Qij = y—{pi — Pj)"^, which 

^ The Coulomb correction applied to a particular triplet is the product of the Coulomb corrections used for the 
three pair combinations contained in that triplet. 
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can be fitted by a double exponential function 



C3 = 1 + Ai exp[-2Q3i?i] + A2 exp[-2Q3i?2] 

witli fitted parameters Ri = 5.08 ± 0.26 fm, Ai = 3.12 ± 0.27, R2 = 1.66 ± 0.08 fm, A2 = 
0.341 ± 0.046 and x^/d.o.f. = 1.10. Such non-Gaussian behaviour was in fact predicted by 



a final-state rescattering model [40|. After the measurement of C2 and C3, the data are 
analysed again and the experimental value of W is calculated using 



W 



{CsiQs) - 1} - {g2(Ql2) - 1} - {C2(Q23) - 1} - {C2(Q3l) - 1} 
2 • V{C2{Ql2) - 1}{C2(Q23) - l}{C2(Q3l) " 1} 



individually for each triplet found, characterized by Q3 and by the values Q12, Q23, and Qai 
corresponding to the three pair combinations contained in the triplet. 
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Figure 17: The factor W as a function of Q3. The error bars include statistical and systematic 
errors. The statistical errors alone are contained within the size of the symbols. 



As described earlier, the estimate of systematic errors is done by varying the different 
analysis cuts, both in the two- and three-pion correlation analysis. The effects on W of 
the statistical errors in the measurement of C2 and C3 are treated as systematic errors by 
changing C2 (C3) by icJCj {^(^Ca)- This last source of error dominates the other ones. All 
these variations are then added in quadrature. 

Figure 17 and Table ^ show VF as a function of foi^ Qs ^ 60 MeV/c, beyond which 
the W significance is too low. The error bars include statistical and systematic errors. The 
statistical errors (not shown separately in Fig. |l^ are by comparison negligible. In view of 
the errors, no significant Q3 dependence is observed. The genuine three-pion correlation is 
found to be substantial with a weighted mean over the five bins {W) = 0.735 it 0.004(stot.) it 
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0.146(syst.). Q 

This result is in agreement with the previously published WA98 result obtained using 
about half of the present data sample |41], where a more detailed description of the three- 
pion analysis method can also be found. More recently, the NA44 experiment [}42[ obtained 
with lower statistics a factor (W) which agrees with our results. 



7 Conclusions 

We have studied mj- distributions for identified vr^ and K~ produced in central Pb+Pb 
collisions at 158^4 GeV. The W.-H. hydro dynamical model has been fitted to the pion spec- 
trum. The resulting fitted parameters favor a combination of a relatively low temperature 
T ~85 MeV and an average transverse flow velocity (/3t) ~0.50. The shape of the vr" 
rriT distribution is in good agreement with the vr*^ mx distribution measured in the same 
experiment. 

Bose-Einstein interferometry of tt~ pairs gives fitted radii of typically 7 fm. This has 
to be compared to the equivalent rms radius of the initial Pb ion of 3.2 fm, indicating an 
expanded emission volume at freeze-out. 

The analysis of two-pion correlations has been performed as a function of /ct using two 
different parameterizations in the LCMS. The results are consistent between the standard 
3-dimensional Pratt-Bertsch fit and the Yano-Koonin-Podgoretskii fit. 

A clear dependence of all radius parameters on is observed, with a stronger dependence 
for the longitudinal radii, indicating a larger longitudinal than transverse expansion. Both 
the cross term Rout-long fro™ the PB fit and vl from the YKP fit deviate from zero, which 
suggests that the source seen within the acceptance does not undergo a strictly boost invariant 
expansion. Moreover the short duration of emission disfavours any long-lived intermediate 
phase. 

A comparison of the data with a hydrodynamical simulation based on the Wiedemann- 
Heinz model, and taking into account acceptance and resolution effects has been made for 
the radii parameters as a function of kx- The agreement is found to be very good. 

The average pion phase-space density at freeze-out has been calculated from measured 
quantities as a function of pT- The results indicate a clear deviation from a Bose-Einstein 
distribution at high pT, but are very well fitted by the Tomasik-Heinz model. The pion 
chemical potential, which is included in the model, is found to be compatible with zero, 
while a large pion chemical potential close to the condensation limit of ttt-tt, seems to be 
excluded. 

Finally, we have studied the 7r~7r^7r^ interference and found a substantial contribution of 
genuine three-pion correlations in central collisions. For <60 MeV/c a weighted mean of 
the strength of the genuine three-pion correlations {W) = 0.735 it 0.004(stat.) it 0.146{syst.) 
was extracted. This is somewhat smaller than what is expected for a fully chaotic and 
symmetric source. 

^ The weighted systematic error is obtained by calculating the weighted average over the five Q3 bins separately 
for each kind of systematic error. These errors are then added in quadrature. On the other hand, adding quadrat- 
ically the systematic errors of the five Q3 bins, as done for weighted statistical errors, would give ±0.078 instead 
of ±0.146 for the systematic uncertainty. 
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Table 1: 3- dimensional analysis as a function of for the PB and the YKP fits. A after correction 
for backgrounds from misidentified pions is also shown. 



{kr) 


0.02 GeV/c 


0.07 GeV/c 


0.125 GeV/c 


0.175 GeV/c 


0.285 GeV/c 


Rts 


6.48±0.18 fm 


6.67±0.15 fm 


5.55±0.24 fm 


5.13±0.32 fm 


4.81±0.29 fm 


Rto 


6.95±0.20 fm 


6.62±0.16 fm 


6.70±0.29 fm 


6.29±0.43 fm 


6.07±0.50 fm 


Rl 


8.69±0.26 fm 


7.54±0.19 fm 


5.95±0.27 fm 


5.45±0.37 fm 


5.15±0.34 fm 


out— long 


1.51±2.11 fm2 


0.55±1.45 fm^ 


-0.52±1.96 fm^ 


-3.01±2.31 fm^ 


-4.51±1.94 fm2 


X 


0.334±0.012 


0.362±0.010 


0.337±0.018 


0.324±0.026 


0.391±0.036 


^cor 


0.646±0.023 


0.634±0.018 


0.557±0.030 


0.511±0.041 


0.596±0.055 


Rt 


6.66±0.14 fm 


6.54±0.13 fm 


5.99±0.19 fm 


5.78±0.29 fm 


5.17±0.26 fm 


R4 


8.68±0.26 fm 


7.52±0.19 fm 


5.97±0.28 fm 


5.49±0.39 fm 


4.94±0.35 fm 


vl 


0.02±0.08 fm 


-0.04±0.06 fm 


0.08±0.08 fm 


0.14±0.10 fm 


0.20±0.09 fm 



Table 2: Average tt phase-space density (/) at freeze-out as a function of pr- 



(Pt) 


0.02 GeV/c 


0.07 GeV/c 


0.125 GeV/c 


0.175 GeV/c 


0.285 GeV/c 


if) 


0.319±0.018 


0.272±0.012 


0.237d=0.019 


0.182±0.022 


0.088±0.012 



Table 3: Weighted mean of the strength of the genuine three-pion correlations (W) as a function 
of ^3- 



Qs 


0.01 - 0.02 GeV/c 


0.02 - 0.03 GeV/c 


0.03 - 0.04 GeV/c 


0.04 - 


0.05 GeV/c 


0.05 - 0.06 GeV/c 


(W) 


830 

-0.175 


780 +°-^^^ 
-0.143 


„ +0.145 
-0.138 


0.654 


+0.177 
-0.193 


557 

"■^^^ -0.293 
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